function [phiDot, torqIncrement] = elevEOM(phi, designVector, setPoint, windProj, contTor)
% elevEOM.m - State equations of motion for the solar panel

l = designVector(3);
w  = designVector(4);

g = 9.81;   %meters/sec^2
t = 0.0005; %meters
rhoAir = 1.2133; %Standard atmosphere Raleigh
u = windProj;
rhoP = 24454.5; % Denstity of monocrystalline silicone plus 5% for framing

%[total_torque, final_velocity, final_Torque] = Controller_SolarDetector(initial_position, final_position, initial_velocity, initial_Torque, cycleTime)

[~, ~, controlTorque] = Controller_SolarDetector(phi(1), setPoint, phi(2), contTor, 10);
torqIncrement = controlTorque;
a = 4*w^2+t^2;
phiDot(1) = phi(2);
phiDot(2) = 6*w*g*cos(phi(1))/a + 3*pi*phi(1)*rhoAir*u^3/(rhoP*t*abs(u)*a)+...
    12*controlTorque/(rhoP*t*l*w*a);